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| Abstract. Linear chains of quantum scatterers are studied in the process of 

lengthening, which is treated and analysed as a discrete dynamical system defined over 
the manifold of scattering matrices. Elementary properties of such dynamics relate the 
transport through the chain to the spectral properties of individual scatterers. For a 
fS |' single-scattering channel case some new light is shed on known transport properties 

of disordered and noisy chains, whereas translationally invariant case can be studied 
analytically in terms of a simple deterministic dynamical map. The many-channel 
^ ■ case was studied numerically by examining the statistical properties of scatterers 

O" 1 that correspond to a certain type of transport of the chain i.e. ballistic or (partially) 

localised. 
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O . 1. Introduction 

J^Tj ■ Linear chains of scatterers represent a useful model of macroscopic structures like multi- 
5_i ■ layered structures, real-life wires, nano-tubes etc. This is the reason that the scattering 



chains have attracted a lot of scientific attention from the early 1980s to the present time. 
Past studies mainly focused on the chain of randomly chosen scatterers and its average 
properties. For good reviews on this topic see [1, 2, 3, 4], where the main approach 
of theoretical analysis is the transfer matrix formalism [5]. Some other interesting 
articles discussing linear chain of scatterers with some disorder using the same approach 
are [6, 7, 8, 9, 10]. There were also attempts to understand the chains using purely 
scattering matrix approach, although they have only partially used its advantages. For 
original references along these lines see Refs. [11, 12, 13]. A considerable breakthrough 
in understanding of localisation in disordered wires was made using the Dorokhov- 
Mello-Pereyra-Kumar (DMPK) scaling equation [14, 15], which gives the scaling of the 
distribution of transmission-like-quantities of individual modes with the chain length. 
More recently there was some renewed interest in scattering formalism, for example in 
the stability analysis of scattering matrix merging procedures [16], and the random walk 
in the scattering chains [17]. 
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In this paper we examine linear chains of abstract quantum scatterers. The 
scatterers on the chain can be any conservative open quantum systems generating 
wave-dynamics, and having two identical waveguides attached on the left/right side 
which connect to the previous/next scatterer (see figure 1). We shall leave out from 
discussion all geometric parameters of the systems and the wave-guides and work 
purely algebraically to make results as general as possible. We are interested in general 
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Figure 1. Schematic picture of a single scatterer representing a basic cell of the chain 
in discussion. 

on-shell transport properties of these chains in the process of lengthening, namely in 
what we shall refer to as the dynamical approach to scattering. We work within the 
scattering matrix formalism in order to avoid divergences present in e.g. transfer matrix 
formalism, and to take advantage of the compactness of the phase space manifold of 
scattering matrices. The presented dynamical approach enables intuitively clear insight 
into the finite and infinite chains. In the following sections, besides discussing dynamical 
properties of lengthening, we also derive some interesting transport properties of the 
chain in three different cases: translationally invariant chain, chain with weak disorder 
and chain with strong disorder. Some of our results can be understood as re-derivation 
of known transport properties of one-dimensional lattices from a simple new dynamical 
perspective, in particular for the case of single- channel scatterers. However, in case of 
multi-channel scatterers we report some new - numerical - results on the scaling of Haar 
measures of scattering matrices corresponding to ballistic and localised dynamics. 

2. Scattering matrix formalism 

In this section we give a short sketch of scattering formalism in quantum chains, but 
the reader can see [18] for details. The stationary Schrodinger equation for the open 
quantum problem composed of the scatterer and two infinite wave-guides is written 
symbolically as 

H\1>) = E\1>) , (1) 

where H represents the Hamiltonian of the system. The wave functions in the left 
and in the right waveguide |t/>r) are expanded in the channel basis {|e^)}^ =1 

d 

I^L,R) = ^^ R |e+) + ^ R |e n -), (2) 

71=1 
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where superscript + and — corresponds to the phase propagation from left-to-right and 
right-to- left, respectively. The number of basis mode functions involved in the expansion, 
denoted by d, is also called the number of (scattering) channels. The smoothness 
of the wave function on the boundary between the wave-guides and the scatterer 
gives the condition which connects the wave functions on both sides of the scatterer. 
The connection between the two sides expressed in vectors of expansion coefficients 
a L,R = {a^' R }^ = i and b L ' R = {& L,R }^=i can be given in the form of a scattering matrix 
S or a transfer matrix T: 

Si>m = ^out , ^ir 
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(4) 



with superscripts L and R corresponding to the left and the right side of the scatterer. 
Following the definitions (3) and (4) it is convenient to express the scattering matrix S 
and the transfer matrix T as block matrices 



S 



r L t R 



T 



X\ x 2 
X 3 Xi 



(5) 



where r L,R represent the reflection matrices and t L,R the transmission matrices for 
incidence of the left and right side of the scatterer. The matrices of this form are 
sometimes called two-way (port) scattering and transfer matrices, because they connect 
two distant openings of the scatterer. Purely from definitions (3) and (4) we obtain 
explicit relations between the scattering matrix and the transfer matrix 



T[S] 



_^R~ 1 r L 



^Ry-R 1 



S[T] 



- x 2 x^x 3 x 2 x~l x 



(6) 



From conservation of probability currents it follows (see e.g. [5]) that the matrices S 
and T fulfill the following relations 



&S = SS 1 = 1 , T^KT = TKT ] = K 



K 



1 









(7) 



implying that the scattering matrix S is unitary, S £ U(2d), and the transfer matrix 
T is hyperbolic, T £ U(d,d). Then by introducing transport probability matrices 
corresponding to the reflection U x and transmission 

Tl x = r^r x , Z x =1?h x , n x + S x = l, x = L,R. (8) 

we define the two basic measures of transport of the scatterer : (i) the average 
transmission (reflection) probability T (TV) as 

TZ=(U X ), T=(L X ) = 1-K, x = L,R. (9) 

and (ii) the standard deviation of the transmission (reflection) probability o\ (cr^) as 



a 



R 



X 



L,R, (10) 
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where we have used the symbol (•) = ^tr{«}. We see that only two quantities can be 
independent and so we use either (1Z, <tr) or (T, <7t) as a measure of transport. 

We build the chain of scatterers recursively by connecting additional scatterers to 
one end of the chain. This procedure is illustrated in figure 2. The chain's length is 
measured in units of the number of scatterers composing the chain. An existing chain 
of length n with the scattering matrix S n is extended with an additional elementary 
scatterer described by a generating scattering matrix S, forming a chain of length n+ 1, 
and with the scattering matrix S n+ \: 



S 



t L 



R 



r 



ti 



71 

R 



r 



The recurrence relation for the scattering matrices of the chain then reads 

Sn+l = S n S , 



(11) 



(12) 



where we introduced a binary operation representing concatenation of the scattering 
matrices. It is explicitly written as 



71+1 



n 



t 



t LR L /-l t R 



n+l 
t R 



t L 



A t L 
i?L ,_1 t R 



(13) 

' n+l ~ 1 ~~<~ L ' n-^ ' °n+l ~ L n- LJ L ) (1^) 

with matrix expressions L = 1 — r^r L and L' — 1 — r L r R . Note that unitary scattering 
matrices U(2d) form a group [19] with the operation 0. We think of the recurrence 
(12) as a discrete dynamical system defined over the space of scattering matrices. We 
study two types of chain-generation: either we take a fixed (static) generating scattering 
matrix S, or S is taken randomly during the growth of the chain. The iteration (12) can 
be written in the transfer matrix formalism in terms of matrix products as 



T 



n+l 



(15) 



where T = T[S] is the generating transfer matrix and T n = T[S n ] is the transfer matrix 
of the chain of length n. But as can be read from the relations (7), matrix elements of 
T n are not bounded in size. This implies that the map (15) is generally unstable and so 
numerically of limited use. 
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Figure 2. Schematic picture of a linear chain of scatterers. It is build by recursive 
lengthening of the initial chain S\ with scatterers S. Scatterer S can be taken as fixed 
or n-dependent. 
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3. Chain of single-channel scatterers 

Here we discuss chains of linearly connected single-channel (d = 1) scatterers. We 
describe the chain and the individual scatterers in terms of 2 x 2 unitary scattering 
matrices S parametrised as 

A = Vl - B 2 , (16) 

where A and B are square roots of reflection and transmission intensity, respectively, 
a L is the reflection phase, and ft ,L are the transmission phases. Therefore any unitary 
2x2 matrix represents a physically legitimate scattering matrix. Note that T[S] = B 2 , 
7Z[S] = A 2 and = = 0. In introduced parametrisation (16) the scattering matrix 
S n of the chain of length n and the generating scattering matrix S read 

S n = S{A n , ofr, ft, ft), S = S(A, a L , ft, ft) . (17) 

The recurrence relation/map for the chain generation (12) can now be written out 
explicitly. We discuss two cases of chain generation. In the first case the chain is 
translationally invariant with fixed S and in the second case S is chosen randomly 
at each iteration step creating a random (disordered) chain. The appropriate explicit 
form of the map differs a bit for the two cases. 



S(A,a L ,ft,ft) 



Ae iah Be il3R 
Be i(3L _ y 4e i(/3L+/3R_aL) 



3.1. Static generating scattering matrix 

Here we consider an initial scatterer with scattering matrix S\ and the generating matrix 
S which is constant along the chain. The infinite chain represents a quantum particle in 
a one-dimensional periodic structure on a half-line with a given initial condition. For a 
scalar periodic potential the problem is discussed in a standard literature on quantum 
mechanics using the transfer matrix approach, see e.g. [20] p. 367. In a related problem 
on a doubly-infinite line, the spectrum has a form of energy bands and the standard 
Bloch theorem applies. In such a case, the explicit form of the recurrence relation for 
the chain generation (12), by using a new variable Xn = ft + ft — (ft + a, reads 



.4 



n+l 



A 2 + A 2 + 2A n A cos x% 



1 + 2A n A cos Xn + (A n A) 2 

Xn+i = Xn + 2A + arg{(l + A n Ae~^){A n + Ae~ [ ^)} , (19) 

ft +1 =ft + ft + arg{(l + A n Ae~^)} , (20) 

ft +1 =ft + ft + arg{(l + A n Ae- i ^)} . (21) 

where we have introduced a parameter A = {ft + ft)/ '2 of the generating scattering 
matrix. The relevant parts of this system of discrete equations are expressions (18) and 
(19) representing an autonomous two-dimensional dynamical system. Note that this 
dynamical system would not change if we restrict ourselves to the case of time-reversal 
invariant (symmetric) S matrices for which ft = ft. Depending on parameters of the 
generating matrix S, namely A and A, we distinguish two topologically different types of 
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dynamics: (i) quasi-periodic or ballistic with an elliptic fixed-point and (ii) convergence 
to a single attractive fixed-point, called localisation. The phase-space portraits for both 
types of dynamics are plotted in figure 3. The study of our two-dimensional dynamical 
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Figure 3. The phase space portraits of the dynamical system (18,19) generated at 
parameters A = 0.5, A = 0.628319 (a) and A = 0.5, A = 0.314159 (b). In (b) we 
connect consequent phase space points in order to give some rough impression on the 
direction of dynamics. Positions of fixed points are indicated by red circular dots. 



system shows that the type of dynamics depends on the sign of the following discriminant 

D = A 2 — sin 2 A , (22) 

yielding for D < the ballistic and for D > the localised dynamics. The expression D 
is a discriminant in the eigenvalue problem of the transfer matrix T[S] (det(T[S]—n 1) = 
0): 



« 1)2 = [cos A ± VD) (23) 

and it is easy to see that, suitably ordering Ki^, 

D < =>• K\ = n* 2 and |ki )2 | = 1 , D > =>• < 1 < |« 2 | ■ (24) 

We can conclude that a growing chain in the transfer matrix formalism (15) for D < 
will be numerically stable and for D > some matrix elements of the transfer matrix 
T n [5 n ] will diverge. 

We note that for the corresponding problem on an infinite line the ballistic case 
D < 0, and the localised case D > 0, correspond to physical energy being positioned 
inside the (Bloch) energy band, or in the gap, respectively. 

In case of ballistic dynamics, D < 0, we find a unique elliptic fixed point (A e , \e)'- 

/—a sin A . , 

A c = u- Vu 2 - 1 , u = , (25) 

71 

Xe = A + - + mvr , (26) 

where m G Z is such that Xe ^ ( — 7r/2, 7r/2). Further, phase-space portrait of ballistic 
dynamics indicates an existence of an additional integral of motion F of dynamics (18, 
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19), which can be derived using the correspondence between the scattering and the 
transfer matrices given by (6): 

P( A ~, \ sinA + A n Asin(A-Xn) , 97 , 

*{An, Xn) = 1 _ A 2 • \ Z7 ) 

By considering the last expression we conclude that the scattering matrix of the chain 
evolves on a one-dimensional manifold. 

In the case of localised dynamics, D > 0, the whole phase-space converges to a 
single fixed point attractor (A a , % a ): 

A a = 1 , (28) 

r / o i sin A / \ 

Xa = A + argjkt + VI - u 2 \ , u = — — . (29) 

J\ 

Hence, the chain converges with increasing length to a state of zero transmission (perfect 
reflection). The convergence is exponential. This can be best seen by locally expanding 
the dynamics around the fixed-point in variables B = yl — A 2 and Xn = Xa + ^Xn and 
so we obtain 

B n , 5Xn~ N n - (30) 

Note that the norms of all matrix elements converge to limiting values, but the phases 
do not, with the exception of ri. Similar behaviour was encountered for multi-channel 
scattering chains. 

3.2. Noisy generating matrix 

Here we discuss disordered chains that are composed of randomly chosen scatterers. 
Such chains of scatterers are known as random chains or random wires. The scattering 
matrix S n of a chain of length n is constructed from its initial state Si by merging with 
generating matrices S that are randomly chosen on each step of lengthening. To simplify 
the discussion we define three disjoint sets of generating scattering matrices: 



M h ={S e U{2) 
Mi ={SeU(2) 
M m = {Se U(2) 



D<0}, (31) 
D > 0} , (32) 
D = 0} . (33) 

We call .Mb the set of ballistic matrices, the set of localised matrices and Ai m 
the set of marginal matrices. In order to measure the volume of introduced sets we 
use a uniquely defined invariant measure over unitary matrices /ir, called the Haar 
measure [21], and normalised so that fi-n(U(2d)) = 1. The Haar measure of the marginal 
matrices /i/f(-M m ) is obviously zero thereby making this set uninteresting for our general 
discussion. The measures of the other two sets are the following: 

/in(-Mb) = 1 - m(Mi) = \ , (34) 

It is important to think about the role of these sets in the construction of the chains. Let 
us construct a chain of length n in which we use m{n) localised generating matrices from 
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A4\. It is evident that in case the ratio m(n)/n in the limit n — > oo is finite then the 
transmission of the chain converges exponentially towards zero. In the opposite case, 
when the chain is constructed mostly of ballistic generating matrices from Aib, the 
transmission can in the worst case decrease linearly with the length of the chain. 

Now we consider chains of single-channel scatterers in which the generating 
scattering matrix explicitly depends on the position in the chain. In the parametrisation 
(17), the dynamics of the chain's scattering matrix in the process of lengthening is 
determined by the system of discrete equations 



B 



B n B 



11+1 



'n 



1 + 2A n Acos((f) n + a L ) + (A n A) 2 

2A + arg{(l + A n Ae~' l{K+ ^ ] )(A n + Ae'' 1 ^ 01 ^)} 



(35) 



(36) 

Pn+i = Pn + arg{(l + A n Ae~^ +a ">)} , (37) 
Pti =Pn+P R + arg{(l + A n Ae-^ +a ^)} , (38) 
where we introduce additional variable n = + (3^ — = \n — & L . This four 
dimensional dynamical system can be reduced to a two dimensional one. It is described 
by the transmission intensity B n , phase <f) n and evolution equations (35) and (36). Each 
iteration step is controlled by three parameters A = A(n), A = X(n) and a L = a L (n). 
This is one parameter more in comparison to the static case of constant generating 
matrix. Just for illustration, we show in figure 4 an example of the chain generated by 
a noisy scattering matrix S chosen deep enough in the ballistic set Aib- We see that 
the presence of a small noise does not strongly deform the trajectory, namely, it does 
not produce localization. It just induces a small variation around unperturbed trajectory 
given by the integral of motion (27) for the average generating scattering matrix (5*). 
Note a similarity with the Kolmogorov-Arnold-Moser type stability in Hamiltonian 
dynamics. 




Figure 4. The phase space portrait in the ballistic regime for parameters Aq = 0.5 
and Ao = 0.628319 without the noise (small points) and with few sample trajectories 
obtained in the presence of noise in the generating matrix S (large points). The latter 
is defined by A = A + e(, A = A + eC, where e = 0.001 and £ S [— 1, 1] is uniformly 
distributed stochastic variable. 



Let us now discuss a general case of random chains with strong disorder. These 
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chains are built using generating matrices from some set A intersecting the set of 
localised matrices A4\ so that ^n(A H .Mi) > 0. We are mainly interested in the 
transmission properties of asymptotically long chains. Intuitively we expect that the 
transmission will converge exponentially towards zero exhibiting exponential localisation. 
In the limit of long chains the transmission is very small and we can replace the exact 
evolution with the following approximation: 

B n+1 =B n f(B(n),4> n + a(n)), (39) 
4> n+l = 4>n + 2A(n) + 2 arg{l + y/l - B(n) 2 e-^ +a ^} , (40) 

/(*, y) = i , (4i) 



2a/1 — x 2 cos(y) 



x 



2\ 



where B(n), a(n) and X(n) are parameters of the generation scattering matrix S G A in 
n-th iteration step. The matrix S is picked randomly at each iteration step and so the 
difference equations (39) and (40) represent a stochastic dynamical system. By taking 
into account the equation (39) the transmission through the chain of length n can be 
written as a product 

n 

B n = B 1 ]Jf(B(k), ( f) k + a(k)). (42) 
fc=i 

In order to understand the scaling of the transmission with the length n we introduce 
the transmission decay rate I n = ^\og(B n ), which reads 

1 " 

/ n = -Vlog[/(£(fc),0 fc + a(fc))], n>l. (43) 
n z — ' 

fc=i 

We would like to express the distribution of /„ over an ensemble of realisations of the 
chain in the limit n>l, defined as 

P n {I) = (5(1 - / n )) stoch . , (44) 

where (. . .) stoch denotes the average over S G A. To obtain this, we need to know the 
dynamics of the variable <p n determined by (40). The distribution of the position of the 
dynamical system <f) n starting at point <pi is defined by 

p n (0;0i) = W-0n)) stoch . • (45) 

It is meaningful to assume that p n has a limiting distribution independent of the initial 
position 0! that is written as 

p(0)= lim^;0i). (46) 

n— >oo 

In some simple stochastic processes we can analytically express p(<fr), but generally this 
is not the case. If p((f>) exists then it is straightforward to show, using the central limit 
theorem [22], that the limiting distribution of I n is a Gaussian distribution 



W~ ( (I -I) 



p n( j ) = A/7^ ex P o 2 ) » n^>l, (47) 
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with the first / and the second moment o~j given by 

7 = (k>g[/(B, + «]) , ^ 2 = (log[/(B, + a] 2 ) - J 2 , (48) 

where (•) is an average over the stochastic variables B and a, and over distributed with 
the density p(4>) (46). The result (47) implies that the transmission decays exponentially 
with the normally distributed rate (47). A similar result was already reported in [7, 11] 
using scaling techniques and transfer matrix formalism, respectively, that are technically 
complicated compared to our derivation. The theoretical predictions are supported by 
numerical studies of which two examples are shown in figure 5. There is a good agreement 
between the measured distribution of I n and the one predicted theoretically. 




-1.8 -1.6 -1.4 -1.2 -1 -0.8 -0.6 -0.4 -2.2 -2 -1.8 -1.6 -1.4 -1.2 

I I 



(a) (b) 

Figure 5. Theoretical and measured distribution of the average decay rate I n at 
n = 100 for two examples of random chains, where A = 7r/10, B G [0, 1] uniformly 
distributed. Variable a L was uniform over [0, 2w] in case (a) and over [0.5, 0.7] in case 
(b). The theoretical curves are given with I = — 1, a\ = 1.4674 in the case (a) and 
with 7 = -1.56325, aj = 1.19556 in the case (b). 



4. Chain of multi-channel scatterers 

We continue our discussion with the general linear chains of <i-channel scatterers where 
the Lie group manifold of unitary U(2d) scattering matrices is (2c?) 2 dimensional [23]. 
Large phase space dimension of lengthening dynamics makes analytical discussion 
very limited, hence our results that we report below are mostly numerical. The 
scattering matrices of the chain are again generated using the general recurrence 
relation (12). We limit ourselves to translationally invariant chains with fixed generating 
scattering matrices. Further, we restrict ourselves to the simplest case of sampling 
the S matrices from the entire unitary group U(2d) with the corresponding Haar 
measure. For considering systems with time-reversal symmetry or some unitary (e.g. 
geometric symmetry), one could follow a similar approach but should first verify that 
the corresponding sub-space or sub-group of S matrices form a semi-group under the 
operation (12). Here our discussion relies on the transfer matrices in the larger 
extent than in the single-channel case. Let us review some algebraic properties of the 
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transfer matrices. The symmetry (7) yields the following relations between eigenvalues 
and eigenvectors of the transfer matrix T 

Tv = kv <^=> (KirfT = — (Kirf , (49) 

K* 

that for two eigenvalues Ki^ with the corresponding right eigenvectors the symmetry 
(7) yields 

- l^lKvx = . (50) 

This means that for every eigenvalue k, with the right eigenvector v, there is a 
corresponding eigenvalue 1 / k*, with the left eigenvector Kv. In case when the eigenvalue 
lies on the unit circle, \n\ = 1, and is non-degenerate, the right eigenvector satisfies the 
relation v*Kv ^ 0. However, for the right eigenvector v corresponding to the eigenvalue 
lying outside the unit circle we have v^ Kv = 0. 

The dynamical system of scattering matrices S n defined by the recurrence (12) has 
all Lyapunov exponents [21] equal to zero and so it is not chaotic. This becomes evident 
from the following discussion. Let us assume that we have a generating scattering matrix 
S and the corresponding transfer matrix T[S] with the spectrum . . . , n 2 d} an d the 
corresponding right eigenvectors Then we can write the lengthening dynamics of the 
chain in the transfer matrix formalism as 

T n = P diag{e' tin }g 1 P- 1 . (51) 

where the matrix P has columns Vi, P = [fjf^i- We distinguish two types of dynamics 
of the scattering matrix description of the chain depending on the spectral properties 
of the transfer matrix: 

(i) If all eigenvalues lie on the unit circle, = 1, V« the dynamics of S n is quasi- 
periodic. This type of dynamics we call ballistic motion and the corresponding generating 
scattering matrices form a set of ballistic unitary matrices Aib- 

(ii) If there is (at least one) eigenvalue outside the unit circle then the transmission 
T[SVi] of the chain decays exponentially towards some plateau value around which it 
oscillates as the chain is lengthened. Following the definition of the transfer matrix (5) 
and the relation with the scattering matrix (6) the average transmission probability can 
be expressed in terms of the lower-right diagonal block X n := (x^) n of the matrix T n as 

7"n = ^tr{(Xt)- 1 X- 1 } . (52) 

The left eigenvectors of T[S] are the rows in the inverse transition matrix P^ 1 = [uf]^ 1 . 
We write the right eigenvector as vj = [a*, G C 2d and the left eigenvector as 
uj = [d,Vi\ e by introducing the upper halves ai,Q G C d and the lower halves 
f3i,rii G C d of eigenvectors. Then according to (51) the dynamics of the block X n reads 

2d 

X n = £ e« n fytf . (53) 

i=l 

We denote by K, — {i : |ac*| > 1} the set of indices corresponding to eigenvalues outside 
the unit circle. The vectors itj, Vi with i G /C satisfy the identity v\Kvi = u\Kui, which 



Dynamical approach to chains of scatterers 



12 



implies that upper and lower halves of these vectors are non-trivial: ||ai|| 2 = ||A|| 2 7^ 
and ||0|| 2 = H^ill 2 7^ 0. We introduce projection matrices P R and P^ onto the set of 
vectors {AjieK an d {Vi}ieic, respectively. Then the transmission can be expressed as a 
sum of non-decaying (oscillating) and decaying term: 

T[S n ] = %[S n ] + 0(exp(-/3In)) , I = logmax{|« m |} , (54) 

rn 

with (3 G {1,2}, where the non-decaying term is written as 

%[S n ] = tr {(Xt)- 1 !" 1 } , X n = (1 - P u L )X n (l - P U R ) , (55) 

where we introduce the decay rate / analogous to that in the single- channel case. If the 
off-diagonal blocks of the matrix X n expressed in terms of introduced projectors are 
zero, (1 - P u L )X n P u R = P u L X n (l - P R ) = 0, then the coefficient (3 = 2, otherwise (3=1. 
The latter is statistically more likely situation. In case the block X n is approximately 
a random matrix then ^[S^] ~ #JC/d, where # denotes the number of elements of a 
finite set. The presented dynamics of transmission is called localised motion and the set 
of corresponding generating scattering matrices are called localised matrices denoted by 
A4\. In case that all eigenvalues are out of the unit circle #/C = d, (P U R,L = 1), the 
transmission decays to zero, and we are talking about total localisation, whereas the 
general case 1 < #K. < d — 1 is referred to as partial localisation. 

The set of generating (unitary) scattering matrices is split into the set of ballistic 
.Mb an d localised matrices M.\. An interesting subset of localised matrices are that 
corresponding to the total localisation. These matrices are named totally localised 
generating scattering matrices and their set is denoted by _Mj* C M.\. The separation 
between different sets of matrices is done on the ground of eigenvalues of the 
corresponding transfer matrices and it is interesting to know the Haar measures of 
these two sets. The measures are obtained numerically by generating unitary matrices 
uniformly with respect to the Haar measure [21] and checking the spectrum of the 
corresponding transfer matrix for eigenvalues outside of unit circle. The result is plotted 
in figure 6. We see that the measure of ballistic matrices decreases very fast with the 
channel number d. We can accurately fit numerical data with an empirical formula 

l^u(M h ) « ttexp(-cud(d+ 1)), with ft w 0.7877, u w 0.4658 . (56) 

In addition we observe that the measure of totally localised matrices also decreases with 
d, however slower, perhaps with a stretch-exponential law. Again we find very accurate 
empirical formula 

I*h(M*) w Zexp(-CVd), with Z w 0.711, C ~ 1-034 . (57) 

From these results we conclude that probability of partial localisation is quickly 
converging to 1 as d increases, however we have at present no theoretical explanation 
or derivation of (56,57). 

It is also instructive to study the eigenvalues of transfer matrices corresponding to 
localised scattering matrices A4\ as they give information about the length-scales of the 
transmission decay T[S n ]. Let us write the spectrum of the transfer matrix T[S] as 

E[S] = {k : det(T[S] - «1) = 0} , (58) 
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Figure 6. The Haar measure of ballistic unitary matrices Uni-Mb C U(2d)) (a), and 
totally localised unitary matrices C U(2d)) as a function of number of channels 

d. The fitted functions arc f h (x) = 0.4658a; - 0.2387 and ft(x) = 1.034a; - 0.342 with 
x = d(d + 1) and x — s/d, respectively. 



and denote by d u (S) the number of eigenvalues in outside the unit circle. We 
investigate the distribution of the maximal eigenvalue modulus: 

P max (t; d) = [ d m {S)S{t - max \E[S}\) , t>l, (59) 

JU(2d) 

and distribution of the relative number of eigenvalues outside the unit circle: 

P u (t;d)= [ dfi u (S)5(t-d u (S)/d), te[0,l], (60) 

JU(2d) 

over the set of scattering matrices S e U (2d) with respect to the Haar measure /in- From 
the distributions -P max and P u we can learn about the decay rates and the percentage 
of unstable eigenvalues involved in the decay of transmission by lengthening of the 
chain, respectively. Both distributions are numerically calculated for several values of 
d and shown in figure 7. In figure 7. a we see that P max (^) has algebraic asymptotics, 
P max (t; d) ~ t~ 3 as t — > oo, and is zero on the unit circle t — 1. By closer inspection one 
can find that the distribution P max (^) has an interesting semi-classical > 1) scaling 
property 

lim VdP max (Vdt; d) = P(t) , t > , (61) 

with P(t) having asymptotic algebraic dependence 

P(t)xar 3 , a = 4.0 ±0.02, t -> oo . (62) 

This scaling of P max (t) implies that the average decay factor increases with the channel 
number d as t ~ \fd and consequently the decay rate (inverse localisation length) 
increases as 7 ~ |log(rf). The scaling law does not work for the single- channel case 
(d = 1), but this is not very surprising. From results for the distribution of relative 
number of unstable directions P u , shown in 7b, we see that the dimension d u on average 
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(a) (b) 

Figure 7. The distributions of maximal eigenvalue modulus of the transfer matrices 
fmax(£ > 1) (a), and the distribution of the relative number x = d u /d of eigenvalues 
outside the unit circle P u (x) plotted against 1 — x (in log scale) (b). The inset in fig. 
(a) shows V / dPmax(v / rfi) as a function of t for several large d, in log-log scale. 



increases with increasing channel number d and the average of the distribution moves 
towards the border value d. We conclude that on average the transmission for larger d 
decays faster and to a lower asymptotic plateau given by TofS^] (55). 

5. Summary and Conclusion 

We have devised an alternative strategy for the analysis of linear chains of scatterers. 
Our approach is based on defining the lengthening of the scattering chain as a dynamical 
system, and connecting its dynamical properties to physical (transport) properties of 
the chain. 

The lengthening dynamics has been shown to be stable for arbitrary number of 
scattering channels. We have been able to reduce the single-channel case to a simple 
two-dimensional dynamical system in which we obtained many results analytically, for 
example the known regimes of ballistic or localised transport correspond to elliptic or 
attractive fixed point of lengthening dynamics, respectively. In the general multi-channel 
case we were only able to give numerical results. We again separate the motion into 
ballistic and localised on the bases of the generating scattering matrices. In particular 
we give accurate numerical results on the scaling of Haar volumes of the sets of scattering 
matrices corresponding to totally localised, and ballistic transport. We have shown that 
in the regime of large number of scattering channels, the transport is most likely to 
be partially localised. In addition, we have examined the distributions of eigenvalues 
of transfer matrices with respect to the Haar measure for the corresponding scattering 
matrices. Interestingly, we found that the distribution of maximal eigenvalue modulus 
satisfies a scaling relation in the regime of large number of channels, corresponding for 
example to semiclassical situations (such as e.g. studied in [24]). 
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